SUNY-NTG-94-35 



General Correlation Functions in the Schwinger Model 
at Zero and Finite Temperature 



James V. Steele, Ajay Subramanian, and Ismail Zahed 

Department of Physics 
SUNY, Stony Brook, New York 11794-3800 



Abstract 

The general correlations between massless fermions are calculated in the Schwinger model 
at arbitrary temperature. The zero temperature calculations on the plane are reviewed 
and clarified. Then the finite temperature fermionic Green's function is computed and 
the results on the torus are compared to those on the plane. It is concluded that a simpler 
way to calculate the finite temperature results is to associate certain terms in the zero 
temperature structure with their finite temperature counterparts. 
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1. Motivation 

A testimony to the complexity of QCD is in the many outstanding problems of the 
theory that are still not solved. For example, it is not fundamentally understood how 
quarks are confined or chiral symmetry is broken, though nature seems to adhere to these 
properties. Many have studied the Schwinger model (QED in two Euclidean dimensions 
with a massless fermion) J2I, |3], Q in the hopes of gaining intuition for tackling various 
such problems. 

Schwinger [|]] proposed this model to show that gauge invariance of a vector field need 
not imply the existence of a massless particle. The gauge field phase shifts the right and 
left handed massless fermions with respect to one another causing a breakdown of chiral 
symmetry. This leads to an anomaly which takes form as the generation of a photon mass 
(m = g/ The vacuum has a non-trivial topology analogous to the case in QCD. This 
theory also has charge confinement because of Gauss' law in one space dimension. 

Since the Schwinger model is exactly solvable, having a form for the general n-point 
correlation function would be useful in many calculations. The result has been approached 
on the plane 0, on the sphere ||, and on the torus 0. Compact surfaces are preferred 
in calculations because the infinities are kept under control. Fermions on the sphere are 
difficult to define and for finite temperature considerations the boundary conditions of 
the torus are more natural. Therefore, work on the torus is desirable. 

Interest in finite temperature calculations in field theories have led to the evaluation 
of two- and four-point correlation functions ||. In this paper, we will provide explicit 
calculations for a six-point correlation function at finite temperature. When combined 
with the previous results, the present result yields the generic form for the n-point corre- 
lation function for the Schwinger model at finite temperature. Exact finite temperature 
calculations of correlation function in field theory are of paramount importance for the 
understanding of finite temperature spectral functions, reflecting on the effects of mat- 
ter on the fundamental degrees of freedom. Real- and Euclidean-time correlators offer 
important information on fundamental issues related to symmetry restoration. 

Our paper is organized as follows: in section 2, we give a brief analysis of the two- 
dimensional Euclidean Dirac operator in the plane. In section 3, we consider n-point 
scalar correlation functions in the two-dimensional Euclidean plane. In section 4, the 
same analysis is extended to the Euclidean two-dimensional torus, and the three-point 
scalar correlator is explicitly worked out. The result is extended to the n-point scalar 
correlator at finite temperature. A brief discussion is given in section 5. Some useful 
algebra can be found in the four appendices. 

2. Preparation 

The Dirac operator may be diagonalized with a complete set of orthonormal eigenfunc- 
tions: ipip n = «7 At (9 At — igA^)ip n = \ n ip n . This equation's Green's function is G(x, y) = 

V _ if \ i s an eigenvalue then —A is also an eigenvalue since {75, ip\ = 0. This 

also implies that the 2x2 matrix, G, is off-diagonal; so the trace of an odd number of them 
vanishes. 

The eigenfunctions of ip with zero eigenvalue are called zero modes. They span the on- 
shell part of the fermionic field and are dominant in the classical limit. In general the gauge 
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field may be written as = —e^ u d u ((p + /) + d^p. The last term may be integrated over 
in the path integral with no consequence. The function / may be chosen to strictly carry 
any singularity of the electric field, called a vortex configuration, and <fi is the quantum 
fluctuation around the vortex. They appear in the action as ^F 2 = -(</> + /)D 2 (0 + /). 
From the index theorem 

k = ^- d 2 x E = ^- d 2 x Savdu,A v = ^- / d 2 x □/, 

LlX J III J ' ZTT J 

the background field, /, must satisfy f(x — > oo) = | ln|x|. Therefore the strength of the 
vortex is related to k, the number of right-handed minus the number of left-handed zero 
modes. The fact that this quantity is not zero is due to the anomaly. 

It can be further shown that the zero modes of a vortex are all of one chirality, 
hence \k\ is the total number of zero modes. This decomposition restricts integration of 
the gauge field on the plane in the path integral to integration over and a sum over all 
possible k. The three gamma matrices (70,71,75) can be chosen to be the Pauli matrices 
from which the identity e^ulvls — ^7 M follows. Finally, solving ipx — for a given flux 
on the plane gives 

with z = x° + ix l , p = 1, . . . , I A; I, and N Pt k is the normalization. 
3. The Plane 

The n-point scalar correlator (hereafter called the n-a correlator) may be separated 
into sums of all possible combinations of the chiral operators: 

- — )*(x) = *P±#(x) (0.1) 

and evaluated using the ideas of Bardacki and Crescimanno ||. A general combination 
of these operators evaluated in a flux sector k will be notated by: 

1=1 

where the e^'s can be either plus or minus. Since the Green's function is off-diagonal, 
terms only contribute for n — k even. The transformation *ff(x) = e -975 '^^^) separates 
the bosonic and fermionic integrals. Due to the chiral anomaly, this transformation has a 
non- trivial Jacobian e ™ 2 / rf2 ^ D (i</>+/) ; which 

can be computed by building up infinitesimal 

transformations on \1/ (see Appendix A). 

For the special case of k = n, only zero modes contribute to C" + . In fact, the 
fermionic integral will vanish unless all the zero modes are taken into account here. A 
detailed calculation given in Appendix B gives 

C n + ... + (x u . . . ,x n ) = M (n) e 2 ^e 4 ^?^ 
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= M {n) e 2a2nK -e g J] 

i>j=i 

The last equality can be proven by setting various z^'s and Zj's equal in the determinant 
and noting this is a zero since two lines are then equal. Here K xy is the bosonic Green's 
function. It can be calculated by introduction of an infrared regulator, fi, which will be 
taken to zero at the end of the calculation. 

K 1 ~ 1 ' 1 



,-,/ - —, ^ 2 ( X -y) = — ( n 1 a - ^V 2 (^ - y) 

U z — m z LJ m z v u — m z □ 



lim— ( 1 - - 1 )<5 2 (a:-y) = -— [K (m\x - y\) + In ^-\x - y\ 

is 1 1 m 
A xx = — In 



Ct. + (*i, • • • , x k ) = M {k) ^(^Y (K -\- 2 ^ = CZ k Axu ■ ■ .,*»). (0.2) 



2g z - fj, 

Since \i does not depend on k, it is the same parameter for all Green's functions. Substi- 
tuting in for all terms, the \zi — Zj\ 2 terms from the zero mode determinant cancel with 
the terms from the Green's functions giving 

m k ( 2 \k(*-l) - 2 E^ 

V 

where Kq V = K (m\x - y\) and M (k) = N x 2 k . . . N k 2 k . 

By taking | for all i and j, the cluster approximation reduces this 

correlator to a form only depending on the condensate 0j, = —me' y /27t 

and so a general form for Mq^ may be found. 

■ 2tt- 

On the plane, the p = k zero mode is non-normalizable. This takes the form of an 
infrared (large distance) divergence. However, this is exactly canceled by the infrared 
(small momentum) divergence of (□ — /i 2 ) -1 as the regulator \i — * as seen in equation 

Calculation of non-minimal correlation functions requires knowledge of the fermionic 
Green's function |J of ip = i(0 — g^$f) in arbitrary flux k. For convenience, we choose 
k > 0. 

0G{x, y) = 5 2 (x - y) - P(x, y) 

f 2i[d z - g(d z f)]G-+{x, y) = 6 2 (x-y)- P{x, y) 
\ 2i[d s + g(d s f)]G+-{x,y)= 8\x-y) 

_ g{f(z,z)-f(z',z')} 

G_ + (x,x') = (0\* L (x)V R (x')\0) = 2m{ -_- /) ; = 

g[f(z',z')-f(z,z)] 

g + 4x,x') = (o\* R (x)* L W\o) = 2m(z _ z>) ; g ++ = 0. 
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Here P(x, y) is the projection operator onto the space of zero modes. The zero modes are 

all right handed and therefore P(x, y) appears in only the G |_ equation. For k < the 

above equations would be exactly the same except that P(x, y) would only appear in the 
G+_ equation. In fact, P(x,y) may be ignored in the calculation of correlation functions 
since the product of zero modes and contractions are antisymmetrized as expounded 
upon in Appendix C. In general, an analogue of Wick's theorem may be used to write out 
all contractions of ^rs and \I/l's possible. The chiral neutral terms (like ^r^l) can be 
replaced by the corresponding Green's function above, and the purely chiral term with 
all the ^/fli's (for k = ±|&|) can be replaced by the zero modes. 

Since G + -(x,y) = G*_ + (y,x), the general non- minimal correlator may be written as 
(for r — s = k > and r + s = n): 



^ {+...+)( — ) ("^1 ' ■ ' ■ ) i "^f +1 j • • • > Xj-jf-g) 



^ n 

2g 2 nK xx V E,,,,/\, ( ._ 2g E e z f(xi) 



X+(zi,Zi) ■■■ X+(^i^i) G + -(x 1 ,x r+1 ) ... G + -(xi,x r+s ) 



e ~ <=i " ' x 

2 



(0.4) 



^r) • ■ ■ X+(^V? G-\ — (x r , X r -(-i) . . . G^ — {x ri Xfjj-s^ 

with Ci = + for 1 < i < r and = — for r + l<i<r + s. With this configuration, 



e * e i 

i>3 



(r — s) 2 — (r + s) k 2 — n 



2 2 
Writing out explicitly all the terms, equation (4) becomes 



/m\n, 2 sk 2 -n 1 iL,,,,^''' 2g E etffpn) TT , 



(2vr; 



e -fl/(a!l) Zie -9f(xi) 



«1 — Zr+1 



e -9f{x T ) z e -gf(x T ) z n,- ig 



k-l p -gf( Xr ) e 3[^r+l)-/(*r)] 
•2r — -2?'+ 1 



e ffl/(*r+*)~/( iC l)] 

2^1 .Zg — |_ S 



It is now easy to see how the / dependence from the determinant cancels the extra term 
from the bosonic integration. Also, as in the minimal correlator above, the \z% — Zj\ 2 terms 
cancel as well. Since the only allowed k are such that (— l) k = (—1)", the remaining terms 
may be written as: 



Cf+... + )(-...-)( X l5 • • • ) x Tl 



me ' 

47T 



l>E,., ./ v •?*'•' 



e »>j 



This can be inserted into the correlators to give 
y-qj(x)tyty(y)) = (^) 2 cosh(2i^o 2/ 
'^^f(x)^(y)^(z) 



4 



= (^) 3 [cosh(2i^) cosh^ 2 ) cosh(2## z ) + 
- sinh(2K% y ) sinh(2K x,z ) sinh(2^)] 

n^)> = e e- 2 fi-x- (0 . 5) 

*=i fe}e{±} 



The last equation, ( p.5|) is the general n-cr correlator. The sum for the sets of ej's runs 
over all possible combinations of pluses and minuses. It is interesting to note that this 
expression can be viewed as the partition function of a two dimensional system of n 
fermions at fixed points with charges ±1 and potential Kq V which at small distances 
reduces to the potential studied by Kosterlitz and Thouless ||. This observation was 
originally due to Samuel |L0 . 



4. The Torus 

We would like to study the Schwinger model at finite temperature. Placing the model 
on a torus [|7| with sides [0, L] x [0, (3\ gives the appropriate geometry for such studies. 
Finally, we shall impose the thermodynamic limit by taking L to infinity and obtain the 
finite temperature result. Taking (3 to infinity should reproduce the results obtained on 
the plane for zero temperature. 

The zero modes are more complicated here than on the plane (Jacobi theta functions) 
yet are normalizable in the finite volume. 

<p 9 {x) = (S h(x) e =»w-)+/(«)] e [ {p - 2 - MA] {kz ^ im . 

with f(x) = 'y{x 1 ) 2 , U(x) = e 2m[hoX ° /f3+h ^ 1/L] , 

z = (x° + ix 1 )/^, and r = L/(3 similar to the notation of Sachs and Wipf |7|]. The 
singularity that □/ had on the plane is spread out over the torus so as not to ruin the 
periodicity. The A integration may be replaced by a sum over k and integration over 
as in Appendix B for the plane. In addition, h Q and hi are the 2ir periodic degrees of 
freedom gained by changing the topology of the space and they must also be integrated 
over: 

J DA, . . . = £<• . ) A = £/ d 2 hD<^ . . . e-i/^ 2 = 



k k 

£ / d 2 hD(/ ) ...e-U'l> a2 <f> 



2-Kk 2 /m 2 V 



The h's are defined modulo 1 corresponding to a full revolution around the two periodic- 
ities of the torus. 

The 3-<t correlator may be broken up into four distinct terms considering the equivalent 
contributions from the positive and negative k sectors. 

(#*(x)m(y)m(z)) = 2[C 3 +++ (x, y, z) + C^x, y, z) + C l ++ _(y, z, x) + C l ++ _{z, x, y)]. 
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The general form of these terms may be written conveniently as (for r — s = k > 0): 

C(+.. .+)(— ) 0^1 ' ' ' ' ' "^ r ' ^-V+lj • • • 5 ^V+s) 

V?i(xi) ... (p k (xi) S+-(x 1 ,X r+1 ) ... SV-^lj^r+s) 



Z-Xdet'^fc) 

(y9l(x r ) ... ipk(x r ) — (,T r ,X r+ i) ... 5*-) — (x r , X r + S ) 

The finite temperature Green's function in the presence of flux k satisfies the equation: 

ipS(x,y) = 5 2 (x,y) - P(x,y). 

Appendix C shows that the solution can be written as just the k = Green's function || 
with an /i-independent term multiplying it to ensure correct boundary conditions. 

S+-(x,y) = S*_ + (y,x) = 

= ^e 9( ^ ) -^ )+/(j/) ~ /(:c)) f/(x)f/ t (i/)e- 2 ^ ^ 1 ~ 22) ^ 1 (0)r fc (2 1 ,2 2 ) 
2irip 



where T k (z 1 ,z 2 ) 



$ z {kzi\ikT)$ A {zi - z 2 + H) 



^{kz 2 \ikr)^{H)d 1 {z 1 - z 2 
These theta functions are related to those from Sachs and Wipf by: 

-a 
-b 



and H = hi — irho. 
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(z, it) = e 2ma(z+b+1 2 Ta ^ 3 (n(z + b + ira)\ir) 



and the tt in the argument and the dependence on r has been suppressed for ease in 
reading. 

Continuing the calculation of the 3-cr correlator, it may be simplified after the 
integration to 



Ql f X y Z \ = V2T\v(i T )\ 8 e Gg 2 K xx +4g 2 (K xy -K xz -K yz ) e -27r/rn 2 V 



X 



X / ^ 2 / le -2^ 2 +47r/io(x 1 + ? / 1 -^)//3 



d 2 {z 2 + H) Ty{z 2: z z ) 

Here K xy is the bosonic Green's function as before on the plane. Its explicit form on the 
torus is (see Appendix D) 



Ag 2 K xy = 4vrA m (x -y)+\n 



r){ir) 



#i(zi - z 2 

After also identifying the finite temperature condensate 



(tytyp = 2 l^( 2T )| 2 c -27r/m 2 V c 2.q 2 ^^ 



o 2lt , 1 l.n 4:71 

2 + —{x 1 -y l ) 2 + ——. 
V m z V 
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the final form for finite volume is 



(0.6) 



with 



\M*i ~ Z2)Mzs)\ 2 J \MH)\ 



d 2 h 



-2Tv{x 1 +y 1 -z 1 -Lh ) 2 /V x 



2 (*i + H) >& 3 {zi)Mzi ~ ^ + H)A{z 2 - z 3 ) 

$ 2 (z 2 + H) ^3(^2)^4(^2 - z 3 + - ^s) 

The full massive bosonic Green's function, 27iA m (x — y) may be rewritten in the limit of 
L — > 00 giving : 



-27rA m (x)= Y, K (mJ(x°-f3ky + (xiy) 



which shows the zero temperature contribution is just KQ(m\x\) as obtained on the plane. 
There is also an extra topological factor, F ++ _, which is not able to be exactly integrated, 
but as L — > 00: 



"&i{z) — > 2e 4 sin7rz ^ 2 (z + H ) — > e 



-ttt(j— ho) p iri(z+hi) 



$ 3 (Z) -> 1 $ 4 (Z + H)^l- e ^(2h -l) e 2ni(z+ hl )_ 

and F ++ _(x, y,z) — > 1. Taking (3 — > 00 gives the same result as on the plane, thereby 
confirming our results. 

Considering ways in which the general result could depend on f} such that the zero 
temperature limit is equation ( p.5|) and the 3-cr correlator is equation ( p.6| ) compel us 
to one conclusion. This is that the finite temperature result can be obtained from the 
zero temperature result by simply replacing the zero temperature condensate and Green's 
function with their finite temperature counterparts. If this is indeed the case, the n-a 
finite temperature correlator can then be written as 



X j 




4?r S eiejA m (xi~Xj) 



In the bosonized form of the Schwinger model |TT] the structure of this result is expected 
T2 \. It is reassuring that the exact finite temperature calculation done here confirms this, 



giving us some confidence in the generalization procedure at finite temperature and infinite 
volume. The generalization to finite temperature and finite volume is less straightforward. 



5. Conclusion 



Two important ways to calculate the n-a correlation functions in the Schwinger model 
were looked at specifically to find a general form for the n-point scalar correlators. The 
plane is easier to work on due to the absence of the cumbersome theta functions, but the 
validity of results from the plane have been questioned in the past. Therefore the result 
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was also derived on the torus and shown to be equivalent to the plane in the large volume 
limit. 

In fact, both methods are very similar. The large distance divergence in the zero modes 
are canceled by the small momentum divergence in the bosonic Green's function. The 
final result for the finite temperature n-a correlator is identical to the zero temperature 
result if the zero temperature condensate and Green's function are replaced by their finite 
temperature counterparts. There is no other spurious temperature dependence as could 
be expected due to the differing topologies in which the two problems were solved. 

The results presented here are directly amenable to the spectral analysis discussed 
by Fayyazuddin et al. Q. They show that in the Schwinger model, the four-, six-, 
n-point scalar correlation functions have a screening length that is determined by the 
meson mass m = e/y/n and independent of the temperature. This is not the case of 
the two-point gauge invariant fermion correlator, where the screening mass was found to 
asymptote nT/2 due to the emergence of two-meron configurations at high temperature 
||15|| . While the present calculations are schematic and far from the real world, they still 
provide nontrivial insights into nonperturbative physics at finite temperature. They are 
certainly of some interest for comparison with finite temperature lattice simulations of 
the model @. 
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Appendix A 

The non-invariance of the fermionic measure under the chiral transformation = 
e~ 915 ^ x ^x{ x ) leads to a non-trivial Jacobian. First noticed by Fujikawa in four dimensions 
TM, we give a derivation in the two dimensional case because there are certain subtleties 



and we have not seen an explicit calculation in the literature. 

Expanding ty(x) and x( x ) m orthonormal eigenmodes of the Dirac operator, ip n (x), 
the above transformation for an infinitesimal change, 6<p(x), in the quantum fluctuation 
field gives 



DtyDty = exp 



d x 



x , 



Yl ^n0)75^n0)l DxDx 



_ e m 2 J (PxSifienvdnAv D^Dx 

for the measure. However, A v = —e up d p (5(f) + f) and so 



>□ 



+ /)■ 



Integrating the infinitesimal transformations, the additional term in the Lagrangian for 
finite is m 2 / <i 2 a;0n(|0 + /). 
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Appendix B 



The n-a correlator may be broken into s + and s_ parts and simplified by the chiral 
transformation discussed in Appendix A to separate the fermionic and bosonic integra- 
tions. 

/n 
i=i 

/n fL 
DxDxD^ J] xPeiXfa) e~ * d2 ^Px+J^f^f) e -\j ***(dWd)* 

K 1=1 

where S eff = J d 2 x[-xiPx + \{<t> + f)^ 2 - m 2 D)(0 + /) + ^-faf 

n 

J = (n 2 - m 2 n)f + 2<? £ e t 5\x - x t ) and f)=$- 91$ f. 



i=i 



Now completing the square in the bosonic integral (with B = D 2 — m 2 D): 



Therefore the bosonic Gaussian cancels the same factor in the denominator. Further 
calculation yields: 

- J d 2 xJB- x J = -J d 2 xfBf + 2g^e i f(x i ) + 2g 2 '£e i e j K XiX . 

i=l i,j 
~y r n n 

= -J d 2 xf{U 2 - m 2 U)f + 2g^ eif{xi) + 2g 2 nK xx + 4g 2 ]T e^K XiXj . 



i=l 



Antisymmetrizing the fermionic functions: 



E(-irn^ ) (^)) = M ( Y ) v^/(- 



1 Z l 
1 ^2 

1 z k 



z l 

z 2 



Jfc-1 



»>j=i 



The non-zero mode fermionic integration leads to a factor: 



det($) 



exp 



^ Trln (^ 



since 



~ 2 

Trln(^r) = -Tr /°° -e^' = m 2 / d 2 x/D/ 

v — S ' Jo da Jo s J 
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where f) a = e a9 ^ f $e a9 ^ f . 

Gathering all terms together, the zero mode saturated correlator becomes 

1=1 i>j 



where \zi — Zj\ = \xi — Xj\ was used. 



Appendix C 

We derive the fermionic Green's function in the presence of a non-trivial flux k on the 
torus. This is required for the non-minimal correlation functions at finite temperature. 
Although other authors have tackled this problem | 16fl , a more explicit form is needed for 
practical calculations. The Green's function is defined by: 

ipS(x, y) = 5 2 (x -y)- P(x, y) 

where P(x,y) is the projection operator onto the space of zero modes. For convenience, 
k will be taken positive. 

Using the trivialization of the U(l) bundle on the torus as in Sachs and Wipf 0, the 
gauge potential is explicitly: 

$ i 27r , 
A = -—x + —h - d X (p 

2tt 

Ai = j^h 1 + d Q 4>. 



Noticing that: 



we can write 



= e 9^^S^)]u(x)i$U\x)e 9 ^ {x)+}{x) \ (0.7) 



S(x,y) = e- 9 ^ {x)+s{x)] U{x)g{x,y)U\y)e~ 9 ^ y)+f{y)] . (0.8) 
This reduces our problem to solving the equation: 

0g(x, y) = 5 2 (x - y) - P(x, y) 

where P(x,y) = e 9 ^ {x)+f{x)] U\x)P(x, y)U(y)e 9 ^ {y)+fiy)] . 

Writing the components of this matrix equation out and letting z = (x° + ix 1 )/ (3 and 
w = (y° + iy 1 )/f3, 

—f3d z g^ + (z — w) = 5 2 {z — w) — P(z,w) 
-(3dzg + ^(z - w) = 5 2 (z-w). 
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Following the ideas of ||, we realize that g + - and g_+ both appear as contractions of left- 
and right-handed fermionic fields in calculations of the non-minimal correlation functions. 
Writing 

k 

9-+(z,w) = g^ + (z,w) - Xn(z)ri*(w) 



in \ 

n=l 



with Vn(z)vl(w) = P{z, w), 

n=l 

-pd z Xn(z) = r} n (z), and - (3d z g„ + (z, w) = 5 2 (z-w). 

The summation above is a summation of the zero modes. In the calculation of non- 
minimal correlation functions, the product of zero modes and contractions needs to be 
antisymmetrized. As shown in section 4, the Green's function appears along with all of 
the zero modes in a determinant. For this reason, the extra summation above will not 
contribute. Therefore only g + _ and g_ + need to be calculated. 

In the presence of a non-zero flux, fields at x l ) and (x°, x 1 + L) are related by the 
gauge transformation: 

A^x , x 1 + L) - A^x , x 1 ) = <9 M a, x 1 + L) = e ia ^(x°, x 1 ) 

2irk n 

where a = — x . (0.9) 

From equations CL8 and p.9| , the boundary conditions on g can be found to be: 

g+4z + l,w) = -e- 2mh0 ~g + 4z,w) g+_(z,w + l) = -e 2mh °~g + _(z,w) 

~g + _(z + ir, w) = e^e-^'+^g+^z, w) (0.10) 
g+-(z,w + ir) = e-* kT e 2 * i{kw+hl) ~g + -{z,w). 
The general form of g + „ is 



1 ^(0) 
2m p v 1 (z — w) 



and f(z,w) must be chosen to complete the boundary conditions in equation p,10| , have 
no poles for z = w, and be 1 at z = w. A solution that satisfies these conditions is 

fl 4 (z -w + H) $ 3 (kz\ikT) 2niho(z _ w) 
) MH) # 3 (hv\ikT) 

A similar calculation can be done for g_ + giving finally: 

~ 1 fMUz -w + H)0 3 (kz\ikr) 2Mz - w) 

g+ ~ { ' ] 2m{3 Mz ~ w)# 4 (H)# 3 (kw\ikr) 

5 (zw) 1 ^Mz-w^M^kr) 2mhM 
y +{ ' 1 2m(3 ^(z -w)$ A {H)d 3 {kz\ikr) 
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Inserting this in equation |0.S| , we obtain the fermionic Green's function. 

The reader might be alarmed that we seem to have managed to invert the Dirac 
operator which is singular in the presence of a non-zero flux. If the solution is not unique, 
it will differ from the one above by a linear superposition of zero modes. This additional 
term does not contribute in the calculation of the correlation function by exactly the same 
arguments presented earlier for going from g_ + to 

Appendix D 

The bosonic Green's function, K xy , may be broken up into a massive and massless 
part: 

Kxy = ^ □ (□-m 2 ) ^ = D-m 2 ~ □ = 

= — {A m (x -y) + —j- - A(x - y)) 

as done in ||. 

The massless part may be calculated by using the eigenmodes that span the allowed 
field (f) on the torus. 

^-y) = (Ahv) = -^ Ux)m 



With <f> n {x) = J-^H^/P+n^/L) and ^ m = Snm 

vv 

Since J <ft = is a necessary condition for a one to one transformation of the fields, the 
value n = (no,ni) = is omitted from the summation. The eigenvalues due to the 
periodicity of the torus are 

/27m \2 /27mi\2 

^ = irr~) + (— ) 

as found in Sachs and Wipf @. 

The sum may be broken up into two terms: one with no ^ and the other with no = 
and Hi 7^ 0. Evaluating the first ni sum and relabeling the remaining summation index 

as n, 

. , 1 „ r cos(27m^) f x i x i , T cos(2t™^) , 
A(x) = — — — coth(7rrn) cosh(27rn— ) — smh(27rn— ) H - — r 

Z7T yj^g 71 P p 7T Ti 

Letting z = (x° + ix 1 )/ (3 and using the identities: 

cos(fcx) cosh(%) = Refcos k(x + iy)] and cos(fcr) smh(ky) = Re[— isinfc(x + iy)], 
the first summation in the braces may be written as: 

r cos27rnz / , , , , e 27rm2 
Re y coth (nrn) - 1 + 



n>0- U 71 
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= Re£[ 



n,r>0 



p 2ninz _i_ p —2irinz 

n 



ln(l - e 



2niz y \ 



-Reln{2ie m2 sin(7r^) JJ(1 - 2e" 2 " Tr cos(2tt2) + e~ 4 " rr )} 



r>0 



since coth(7rrn) - 1 = ^ Q e 2nTnr and ln(l — x) — ~^ ^- 
The second summation in the braces gives 



TXT , TT(X 

— nx H 

6 r 



1\2 



The first two terms can be combined with the first summation to form the product 
representation of ~&\ and 7](ir) to finally give 



— An A(x — y) — In 



T]{IT) 



#i(t(zi - ^2)|^) 



2 27T. ! ,. 2 

+ y(^ 1 -y 1 ) 2 . 
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